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The nature of cosmic time is illuminated using Hamilton- Jacobi theory for general relativity. For 
problems of interest to cosmology, one may solve for the phase of the wavefunctional by using a 
line integral in superspace. Each contour of integration corresponds to a particular choice of time 
hypersurface, and each yields the same answer. In this way, one can construct a covariant formalism 
where all time hypersurfaces are treated on an equal footing. Using the method of characteristics, 
explicit solutions for an inflationary epoch with several scalar fields are given. The theoretical 
predictions of double inflation are compared with recent galaxy data and large angle microwave 
background anisotropies. 
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One obtains a better appreciation for the nature of cosmic time by studying the role of inhomogeneities in general 
relativity 0]- Q. After all, a time-hypersurface represents the arbitrary manner in which one slices a 4-geometry. 
y—t , However, there are numerous problems associated with quantum aspects of time 0. One should perhaps be content 
' to consider the semiclassical limit of the quantum theory where one retains only the leading order terms in Planck's 
^sO \ constant Ti that appear in the Wheeler-DeWitt equation. This approximation leads directly to the Hamilton- Jacobi 
(HJ) equation for general relativity. In the present work, superspace solution techniques of the HJ equation will be 
further advanced to include the case of several scalar fields interacting in a cosmological setting. These methods 
simplify the comparison of theoretical models with astronomical observations. 

The problem of choosing a gauge in general relativity has complicated the interpretation as well as the application 
, of the theory For example, when solving Einstein's equations, one must make arbitrary choices for the space-time 
coordinates. Many of these decisions may be postponed or even avoided by solving the HJ equation for general 
relativity. Hamilton-Jacobi theory leads to a covariant description of the gravitational field (l|- ||. It yields a 
I generalization of earlier work on covariant cosmological perturbations by Ellis et al J8| . 

Already, HJ techniques have been applied to numerous problems in cosmology including: (1) a detailed computation 
■ of microwave background fluctuations and galaxy-galaxy correlations [ IL generated in the power-law inflation model 
p| which arises either in induced gravity jl0| or extended inflation |L2f; (2) a relativistic approach to the Zel'dovich 
approximation |l^] describing the formation of sheet-like structures during the matter-dominated era |l4|] , |l5| , [|l6| ; 
(3) attempts to recover the inflaton potential from cosmological observations ju|, ]l7j ; (4) a construction of inflation 
models that yield non-Gaussian primordial fluctuations |j| — - such models may alleviate the problems of large scale 
j_j ■ structure in the Universe |l8| . 

In Sec. II, I set forth the HJ equation and the momentum constraint equation describing the interaction of two scalar 
fields with the gravitational field. The extension to more fields is straightforward. The object of chief importance is 
the generating functional S which is the phase of the wavefunctional in the semiclassical approximation. In Sec. Ill, a 
spatial gradient expansion is used to solve the HJ equation. The zeroth order term for the generating functional, S^ ' , 
contains no spatial gradients. It describes the evolution of long-wavelength fields. The second term, S^ 2 \ contains two 
spatial derivatives. Its evolution equation is a linear partial differential equation which may be simplified using the 
method of characteristics. In particular, by invoking a local change of variables, one may solve this equation by using 
a line integral in superspace. Different contours of integration correspond to different choices of the time-hypersurface, 
but they all lead to the same answer for S^ 2 \ In Sec. IV, it is shown how to effectively sum an infinite subset of 
terms in the spatial gradient expansion. This generalizes the quadratic curvature approximation that was employed 
by Salopek and Stewart |jj . The resulting Riccati equations may be reduced to ordinary linear differential equations 
using matrix methods. In Sec. V, initial conditions are given for a period of cosmological inflation. In Sec. VI, some 
examples and solutions are given for inflation with two scalar fields. Using analytic methods, one is able to compute 
fluctuation spectra that give rise to galaxy clustering. Previously, one required a tedious numerical computation [flo| . 
Theoretical models are constrained by large angle microwave background anisotropies and galaxy clustering data. 

(Units are chosen so that c = 8irG — 8tt /m 2 p = % = 1. The sign conventions of Misner, Thorne and Wheeler |n]] 
will be adopted throughout.) 
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II. THE HAMILTON- J ACOBI EQUATION FOR GENERAL RELATIVITY 



For two scalar fields <j> a , a = 1, 2, interacting through Einstein gravity, the Hamilton- Jacobi equation, 

ss ss 



= H{x)= 7 1/2 [27ik(x)7j/(x) - ~fij( x hki(x)} 
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5(f> a (x) 



a=l Z 

and the momentum constraint equation, 

= Hi(x) =-2 ( 7lfe - 
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(2.1b) 



govern the evolution of the generating functional, S = S[4> a (x) , 7^ (x)] , in superspace. For each field configuration <j) a (x) 
on a space- like hypersurface with 3-geometry described by the 3- metric 7^ (x) , the generating functional associates a 
complex number. Here R denotes the Ricci scalar of the 3-metric, and V[<fi a (x)} is the scalar field potential. In the 
ADM formalism, the space-time line element is written as 



ds 2 ee g^dx^dx" = (-N 2 + j ij NiNj) dt 2 + 2N t dtdx l + j tj dx l dx j 



(2.2) 



where N and N{ are the lapse and shift functions, respectively. The generating functional is the 'phase' of the 
wavefunctional in the semiclassical approximation, 



(2.3a) 



(The prefactor is neglected here although it has important implications for quantum cosmology pl|.) The probability 
functional 



V= |*| 2 , 



(2.3b) 



is just the square of the w avefunctional (see, e.g., ref. [g0|). 

The HJ equation ( 2.1a ) and the momentum constraint ( 2.1b ) follow, respectively, from the Gq and G° Einstein 
equations with the canonical momenta replaced by functional derivatives of S: 



6S 



84> a {x) 



and tt ij (x) 



SS 



(2.4) 



The momentum constraint demands that the generating functional be invariant under an arbitrary change of the 
spatial coordinates |2^] (see also ref. [^9|, p. 1185). The Hamilton- Jacobi equation for general relativity is analogous 
to the Tomonaga-Schwinger equation that yielded a covariant formulation of quantum electrodynamics p3|] . If the 
generating functional is real, the evolution of the 3-metric for one particular universe is given by 



SS 

(ffij - - N jH ) /N = 27- 1 / 2 (2j jk j a - Hjjki) ^— : 



whereas the evolution equation for the scalar field is 



(<j> a ~N^ t )/N = 1 - 1/2 ^ r - 



(2.5a) 



(2.5b) 



Here | denotes a cov arian t derivative with respect to the 3-metric 7^ . Since the lapse and shift function appear neither 
in the HJ equation (2.1a) nor in the momentum constraint (2.1b), these two equations are covariant in that they are 
valid for all choices of the space-time coordinates. Ail gauge-dependent quantities which are associated with the lapse 
and shift appear only in the evolution equations (2.5a,b) for the 3-metric and scalar field. 
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III. SPATIAL GRADIENT EXPANSION OF THE GENERATING FUNCTIONAL 



Superspace describes an ensemble of evolving universes, and its complexity is overwhelming. One may reduce this 
complexity by expanding the generating functional 

oc 
71=0 

in a series of terms according to the number of spatial gradients that they contain. The great simplification that 
arises is that one encounters linear partial differential equations at each order for n > 1. A careful exposition is given 
of the generating functional of order two. This example illustrates many of the essential aspects required in solving 
more challenging problems. The consistency of higher order solutions is also demonstrated. 



A. Generating Functional of Order Zero 

To order zero, the Hamilton- Jacobi equation is 



1/2,0 s ^ (0) SS l ° ] X^l i/ 2 /«S (0) V 



+ l 1/2 V{4> a ) =0. 



(3.2) 



In the full HJ eq.( [2.1a ), the spatial derivative terms appearing in the square brackets have been neglected. They will 
be recovered at the next order. One attempts a solution of the form, 

5<°>fa,(a0,7<i(a0] = -2 / ^x^ 2 H (4> a (x)) . 



(3.3a) 

for the zeroth order generating functional. The numerical factor —2 is chosen so that the function of the scalar fields, 
H = H(cj> a ), corresponds to the usual Hubble parameter in the long- wavelength approximation. Because the integral is 
over cPxj 1 / 2 , this functional is invariant under spatial coordinate transformations, and consequently the momentum 
constraint eq.( 2.1V| ) is satisfied. The HJ eq.(3.2) of order zero holds provided that H satisfies the following nonlinear 
partial differential equation: 



H- 
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(3.3b) 



Since the metric does not appear in this equation, it is often referred to as the separated Hamilton- Jacobi equation 
(SHJE) (of order zero). 



B. Generating Functional of Order Two 

The second order HJ equation is 

It is a linear partial differential equation of the inhomogeneous type. It may be solved using a two step method. (1) 
One firstly applies the method of characteristics to choose an integration parameter or time parameter; the remainder 
of the variables are then frozen in time. (2) Secondly, following the work of Salopek and Stewart (J], one makes an 
ansatz for the second order generating functional which contains all terms with two spatial derivatives. If one wishes 
to be more direct, one may employ a line integral in superspace [|| which elegantly illuminates the nature of cosmic 
time. 

The characteristic equations 
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ld ' n 4?, (3.5a) 



N dt " d(t> a ' 

2H Jl3 (3.5b) 
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are determined by the coefficierits of the functional derivatives, 6S/54> a and SS/S^ij, appearing in eq.( |3.4| ). Here the 
time parameter t as well as the lapse N are arbitrary, but for the sake of simplicity one can safely assume that they 
are local functions of the scalar fields: t = t((j) a ), N = N(<fi a ). The above equations may be simplified by performing 
a conformal transformation, 

7«(a0 =n 2 (cf> a )f ij (x), (3.6) 

where the conformal factor, f2 = f2(</> a ), is a function of the scalar fields, and the conformal 3-metric, , is independent 
of time t. The reduced characteristic equations are then: 

1 d<t> a dH 

= -2— , (3.7a) 



N dt 
1 dlnft 
N dt 



= H . (3.7b) 



There are numerous ways to solve these equations. If one is fortunate to obtain a solution for the Hubble function 

H = H ((f> a , 4> a ) which depends on two homogeneous and time-independent parameters <f> a , a — 1,2, then one may 
integrate these equations immediately: 

J1 ^") = (^r) ' (3 - 8a) 

\d<pi/ 

tlS dH.dH ,„ nl , 

e(<p a ) = ^/^, (3.8b) 

d<j>i d(j> 2 

fa = r»- 2 (0 a ) nj ■ (3.8c) 



Both e(<p a ) and are independent of time: 



~=0, (3.9a) 
= . (3.9b) 



1 dfjj 
N dt 



By differentiating the SH JE (3.31: ) with respect to the parameters, one may verify that these fields satisfy the reduced 
characteristic equations ( (3.7b| ). Note that this solution refers neither to the time parameter t nor to the lapse N — 
they are redundant variables. 

The integration of the characteristic equations may be viewed as a transformation of the fields 

foi,fc,7«)->(*Wi J -). (3.10) 

Utilizing the conformal 3-metric /y instead the original 3-metric 7^ is analogous to using comoving coordinates rather 
than physical coordinates in cosmological systems. Q,(x) = Q,[<p a (x)] can be interpreted as the integration parameter 
or 'time' parameter, whereas e(x), fij(x) are independent of time. The integration parameter could have been chosen 
to be any function of the scalar fiel ds, bu t the choice of t = = f2(</> ) will simplify the subsequent integration. From 
the second characteristic equation ( [3.7b|) , dVl/dt = 1, implies that the lapse is given by N = 1/(QH), and hence the 
scalar field evolves in f2 according to 

9 JA =^ d JL on) 

dn J e QH d(j> a ' v ' ' 

In some instances, it may be difficult to de termi ne a solution for the Hubble function which depends on two param- 



eters. Alternatively, one may integrate eq.( 3.11 ) where e now refers to an arbitrary and inhomogeneous constant of 
integration. 

Using the new variables, the HJ equation of order two may be rewritten as 
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= I~V2 

2 r 



(3.12) 



The method of characteristics has simplified the left hand side considerably. It now remains to express the right hand 
side in terms of the new variables. For example, a conformal transformation of the Ricci curvature leads to 



fl(fi 2 /«) = rr 2 i? - 8rr 5/2 (^ 1/2 ):* 



(3.13) 



where R is the Ricci scalar of the conformal 3-metric fij , and a semi-colon ; denotes a covariant derivative with respect 
to the conformal 3-metric. As a result, one finds the following final form for the HJ equation of order two, 
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= / l/2 



-fill + 1(9, e) n.ifl' 1 + 2m(Q, e) 9,^ + n(fi, e)e, l e l 

2H H9 ' 



(3.14a) 



where l,m,n are functions of (f2, e): 

l(Q,e) 
m(9, e) 
n(n, e) 
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dcf>2 \ ( dcj) : 
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(3.14b) 
(3.14c) 
(3.14d) 



1. Ansatz for Second Order Solution 



An Ansatz, 



d 3 xf/ 2 



j(9, e)R + ku(n, e)VL-iVt A + 2k 12 (9, e)9- i & i + fc 2 2(^, e)e ;i e ; 



(3.15a) 



will be used to integrate the HJ equation of order two. j, k%i, fci2 and k 22 , are functions of (9, e) which are determined 
to be: 



fcn (O, e) 

k 22 (9,e) 
fci 2 (n,e) = 0, 



dQ! 



1 



2H(Q',e) 



io(e) , 



1 



dtt' n(fl',e)+k (e) 



(3.15b) 

(3.15c) 

(3.15d) 
(3.15e) 



where n(9, e) was defined in eq.( [3.14d ), and jo = Jo(e) and ko = ko(&) are arbitrary functions of e. 

I n orde r to construct the above solution, one computes the functional derivative with respect to 9(x) of the Ansatz 
eq.(3.15a): 



-1/2^ 
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de 



( dk 22 2 9kl2 



V 89 
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e-ie ;t - 2kn9:l - 2k 12 e'.l , 



(3.16) 



and matches terms with HJ equation (3.14a) of order two. Since no e'.\ term appears in eq.(3.14a), one finds immedi- 
ately that the cross term k\ 2 vanishes. This fortunate circumstance arose because we chose 9 to be our integration 
parameter. The remaining equations 
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dfl 
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HQ ' 
-I, 



(3.17b) 
(3.17c) 

(3.17d) 

(3.17e) 



lead to the solution given above. Initially, it appears that this system of equations may be overdetermined since fcn 
appears in the three equations ( 3.17b| -3.17d). It is a minor miracle that given the first equation, the latter two are 

automatically satisfied. 

The integration of the HJ equation of order two is a nontrivial result. The ansatz ( [3.15a ) and the prescription for 
determining its free functions represent a complete and explicit integration of the problem. Previously, Salopek and 
Stewart Q had solved the single field system. The new ingredients required for resolving the two scalar field problem 
were: (1) the meth od of characteristics, and (2) choosing f2(x) = f2[</> a (a;)] as the integration parameter. Even then, 
the final form, eqs.(3.15a-e), of the generating functional of order two is valid for all choices of the time hypersurface. 



2. Line Integrals and Potential Theory 



The HJ eq.( |3.14a ) of order two has the form of an infinite dimensional gradient. It may be integrated directly using 



a line integral in superspace. The required aspects of potential theory are first briefly reviewed. 

The fundamental problem in potential theory is: given a force field g l {yk) which is a function of n variables yk, 
what is the potential $ = $(j/fc) (if it exists) whose gradient returns the force field: 

dy~r gl[Vk) 1 (3 - 18) 

Not all force fields are derivable from a potential. Provided that the force field satisfies the integrability relation, 

(3.19) 



Q _ dg l dgi 
dyj dy t 



d d_ 

dyj ' % 



(i.e., it is curl- free), one may find a solution which is conveniently expressed using a line integral 

/ E^Vfo) ■ (3-20) 



c J 



If the two endpoints are fixed, all contours return the same answer. In practice, we will employ the simplest contour 
that one can imagine: a line connecting the origin to the observation point y^. Using s, < s < 1, to parameterize 
the contour, y x = syi, the line integral may be rewritten as 

*(^) = E / 9 j (Vi)+$o , (3.21) 

3=1 Jo 



where $o is independent oiyu- 



3. Line Integral Method for Second Order Generating Functional 
By identifying, j — > x, — * J d 3 x, yj — > Q(x), dyj — dsCl(x), one may integrate the infinite dimensional gradient 



eq. 



CHI ), 



G 



(3.22) 



(2) 

with = sCl(x). One recovers the result fl3.15a| ) guessed earlier provided the 'constant functional', Sq , is given 
by 



4 2) [e(»), /«(*)]= / d 3 xf^ 2 [io(e) J R+fc (e)/«e i e,. 



(3.23) 



which is independent of Because of an integrability condition discussed in the next section, this result does not 

depend on the contour of integration in superspace. 



4- Integrability Condition for Generating Functional 



A line integral proves useful for computing higher order terms in the spatial gradient expansion. 

For n > 1, one finds that S^ 2 ™) satisfies the following inhomogeneous, linear partial differential equation: 



00) S (2n) = ~Tl {2n) {x). 



The differential operator 0(x) is given by 



O x = V-2 



2H lrj 



(3.24a) 



(3.24b) 



The remainder term TZ^ 2n ^(3 



n-l 2 ro(2p) Jc(2n-2p) 

+^-i/2 v V - — — - + v< 2 ™) 



(3.24c) 



is independent of S^ 2 ™) but it contains contributions from all previous orders — it is assumed to be known. The 
superspace potential y( 2 ™) is defined to be 



V™(x) = -\^R + E W /2 l ij Ma tj for n = 1 , 

a=l 

— otherwise. 

By manipulating eq.( |3.24a -d), one may compute the following commutator, 
0{x), 0(y)\ S (2n) = 0(y) Tl {2n) {x) - 0{x) TZ^ 2n) (y) 

= tft( x ) nf n - 2 \x) + 7 <V) nf n - 2 \x')] ^-f{x - x') , 



(3.24d) 



(3.25) 



which assumes by induction that S^- 2 \S^ 4 \ . . . ,S^ 2n 2 ) satisfy eq.(3.24a). The 'integrability condition' of potential 

j 



theory, eq.(3.19), demands that the commutator (3.25) vanish. In the above expression, TC? n 2 ' is the momentum 



constraint evaluated using the generating functional of order (2n — 2): 

&S(2«-2) 



5-fki{x) 



5S {2n~2) 2 6s {2n-2) 

07*1(1) fr[ Ha{x) 



(3.26) 
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We conclude that 6>( 2 ") is indeed integrable provided the term of previous order, <S( 2 ™ -2 ), is invariant under 
reparametrizations of the spatial coordinates: Hj 2n ^ = 0. In general, the integrability condition for the Hamilton- 
Jacobi equation follows from the Poisson brackets Q between the energy densities evaluated at the two spatial points 
x and x': 



{H(x),H(x')} = [Y^x)n j (x)+^(x')H J (x')} ^Ld 3 (x-x') 



(3.27) 



In practice, one would solve the HJ eq. (|3.24a ) of order 2n by changing variables according to eq.( 3.8a -c), and then 
applying a line integral analogous to eq.( [3.22 ): 



S^[fl(x),e(x)J ij (x)} 



d A x 



1 ds 



o a H[fl(x),e(x)] 



TZ^[Q(x),e(x)J t3 (x)] 



+5q 2 "^ [e(x), fij(x)] , with Q(x) — sO(x) . 



(3.28) 
(3.29) 



Once again, the 'constant functional,' 5q 2 ™' [e(x), (x)j depends only on the fields e(x), (x), which are independent 
of Q(x). The constant functional contains all such terms with 2n spatial derivatives which are invariant under 
reparametrizations of the spatial coordinates. 



C. The Nature of Cosmic Time 



For the semiclassical limit, a line integral in superspace goes a long way towards illuminating the nature of cosmic 
time. Different contours of integration analogous to eq.( p.29| ) correspond to different intermediate choices for the 
integration parameter or time parameter. However, they all yield the same answer for the generating functional 
provided that spatial gauge-invariance is maintained at each order of the computation. 



IV. QUADRATIC CONSTANT APPROXIMATION FOR GRAVITY WITH MULTIPLE SCALAR FIELDS 

In the last section, a prescription was given for computing the generating functional to arbitrary order in the spatial 
gradient expansion. For a single scalar field or a single dust field, this has been done explicitly to 4th order or 6th 
order, respectively, in spatial gradients [^). The case for a single dust field is of considerable practical importance 
because it leads to the Zel'dovich approximation for general relativity Q], Jl5| . 

However, for some cosmological applications, a finite number of terms is insufficient. In this section, it will be 
shown how to effectively sum an infinite subset of terms in the spatial gradient expansion. The 'quadratic constant 
approximation' generalizes the 'quadratic curvature approximation' developed in an earlier paper M. 



A. Factoring Out the Long- Wavelength Background 



Beginning with the full HJ eq.( 2.1a ), one first subtracts the long- wavelength background from the generating 
functional: 



S = + T , 5<°> = -2 J d 3 x 7 1/2 H(<t>a) 



(4.1) 



where H satisfies the separated Hamilton- Jacobi equation ( 3.3b ). The functional for fluctuations, !F, now satisfies 

2 

E„9H 5 J- „ TT 5T i in .„ . . , . . . ... 5T 5 J 7 
" 2 ^TZr + 2F ^^— + ^ [2 lu {x)j jk {x) - 7«(a07fcl(a0 x , w T\ 



a=l 



2 



(4.2) 



The first line of eq.(4.2) may simplified if one introduces the change of variables, (0 aJ 7ij) ~~ ¥ (^,e,/y), described by 
eqs.(3.8a-c). Functional derivatives with respect to the fields transform according to 
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(Henceforth, one suppresses the symbols |n, | e , |/ y which denote the variables that are held constant during functional 
differentiation.) 

The HJ equation reduces to 



sr r 1/2 

sn + n 4 H 



[Zfafjk - fijfkl_ 



ST ST 

Shi Sfk 



E 



f-1/2 

2fl 4 H 



d<t> a \sn~ n il3 5f l3 



de ST 
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i> a Se 



5S™ 

sn 



(4.4) 



Tij ojki a=1 

where 5( 2 ) was given in eq.(3.15a). The momentum constraint maintains the same form as before but it is now 
expressed in terms of the new variables (SI, e, /y): 



ST 



' $fkj 



ST 



ST 



ST. 



Hi(x) = -2 [ f ik — ) +_/ Wii +_ ei + _n i = o . 



Sfi 



Se 



sn 



(4.5) 



B. Integral form of HJ equation 



Following ref. g, an integral form of the HJ eq.(4.4) may be constructed using a line integral in superspace, 



T[n{x) ) e(x),f ij {x)}+ / d 3 x 



1 dsf- 1 ' 2 



ST ST 



— ;; [2/ii/jfc fijfkl] . ,. 



d 3 x 



1 ds 2 / -l/2 







where once again Cl(x) = sil(x). 



d</> a \sTl Jl3 Sf t] 



de ST 



b a Se 



= SW[n(x),e(x)J ij (x)} 



(4.6) 



C. Ansatz 



One makes an Ansatz for the fluctuation functional (4.1) of the form, 

T = S {2) + Q, 

where the quadratic functional, 



Q = J d 3 xf^ 



RS n R+ 2RS 12 (D 2 e) + (D 2 e) S 22 (D 2 e) + R i] T R tj --RTR 



(4.7) 



(4.8) 



consists of all quadratic combinations of the fields, D 2 e, R, Rij (but not D 2 £l), which maintain spatial gauge invari- 
ance. D 2 is the Laplacian operator with respect to the conformal 3-metric, e.g., 



D 2 R = D' l D ;i R = R l . tl = r 1 ' 2 (/ 2 fit.. 



One interprets the operator T for tensor perturbations to be a Taylor series of the form 

f(fi, e ,5 2 ) = ^T n (tt, e )(5 2 y\ 



(4.9) 



(4.10) 



and similarly for the scalar operator S a b, which is a two- by- two symmetric matrix. (The full Riemann tensor Riju does 
not appear in the Ansatz because for three spatial dimensions because it may be written in terms of the Ricci tensor 
p5|.) Although the first and fifth terms in eq.(|4.8|) may be combined into a single one, the present form simplifies the 



final evolution equations for S and T (see eqs.(116b), (4.17b)). The Ansatz (4.8) generalizes the quadratic curvature 
approximation considered by Salopek and Stewart Jl[]. 

The perturbation rules are the same as those developed in ref. pi. By first order, we refer to terms such as R, 
D 2 fl, D 2 e, or D 2 R, D 4 f2, D A e, which vanish if the fields are homogeneous; they may contain any number of spatial 
derivatives. Quadratic terms are a product of two linear terms. 

In computing the various functional derivatives, it is useful to note that for a small variation of the conformal 
3-metric Sfij the corresponding change in the Ricci tensor is 



(4.11) 



In the integral form of the HJ equ ation , integration by parts is permitted which simplies the analysis considerably. 
(However, in its differential form ( 2.1a ) one cannot simply discard total spatial derivatives.) In addition, all cubic 
terms are neglected. 



To linear order, the functional derivatives of , eq.( 3.15a ), are 

(2) dj 



R - 2/cn D Z Q 



.f' //" • f iJ ./:;/" 



Se de 

and those of the quadratic functional, eq.(|4.8|), are 



(4.12a) 
(4.12b) 

(4.12c) 



-1/2 



f 



SQ 
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SQ 



0. 



T 



R™ - D 2 Rf ij 
1 







+ 25 12 


(D 2 ey J -D A ef lJ 



+ -D 2 RP 3 - D 2 R l] 
4 4 



/- 1/2 S _Q = 2D 2 
oe 



S 12 R + S 22 D z e 



(4.13a) 
(4.13b) 

(4.13c) 
(4.13d) 



At the present level of approximation, D 2 commutes with any function of (Q,E), e.g., D 2 Sn = SnD 2 , etc. After a 
straightforward computation (which is very similar to that of ref. 0), one obtains the tensor equation 



— — 1 — — lj + TD 2 ) = , (tensor) 



(4.14) 



and the scalar equation 



f)C „ _ ^ i ^ 

— — + 2D A SAS + D 2 C T S + D 2 SC + -B = , (scalar) 
a In i I 2 



(4.15a) 



where A, B and C are operator matrices. (From now on, the operator symbol^will be suppressed.) They are given 
by 



.4 
B 
C 



1 


-- -T 




i 






i 


—>rp 

a 1 b 1 


n 3 H 



Si 



(4.15b) 
(4.15c) 
(4.15d) 



Here T denotes the transpose of a matrix or a vector, and the vectors a±, a 2 and b\, b 2 are 
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, d In ft de 



(4.15e) 



a 2 



, d In ft de 



(4.15f) 



917 \ / dj j \ dj de d In ft dj ^ de 
d<j>i ) \dfl ft/ dedcf>i' d(f>x de 22 d<j> i 



(4.15g) 



»7 = 



<9ft 



<9j de ^dlnftdj <9e 
de d(f>2 ' 902 de 22 902 



(4.15h) 



The tensor equation is a single nonlinear ordinary differential equation of the Riccati type. The scalar equation is 
also of the Riccati type, but the dependent variable S a b is a 2-by-2 symmetric matrix. An additional complications 
arise in the latter case because the coefficients A, B, C, appearing in eq.( [f.l5a ) need not be commuting. 



D. Reduction of Riccati Equations 



The Riccati equations, ( 4.14 ) and ( 4.15a ), may be reduced to linear ordinary differential equations. For the tensor 
perturbation, one defines the Riccati transformation, T(ft, e, D 2 ) — > y(fl, e, D 2 ), 



2D 4 y d In ft 

which leads to the linear equation for a massless field : 



Hrt 3 1 dy ~_ 2 , 
T = — ^— — — — — — D j , (tensor perturbation) 







d 2 y 



d(lnQ) 2 



1 dH 
H d In n 



1 



dy 



dlnQ H 2 VL 



;D 2 y . (tensor perturbation) 



(4.16a) 



(4.16b) 



Because non-commuting matrices appear in the scalar case, one must be careful about the order of various ma- 
trices, but otherwise the method is the same as the previous case. By defining the matrix Riccati transformation, 
S^e,^ 2 )^^^,^ 2 ), 



1 A-i£^W-i 



1 



2D 4 dlnfl 2D 2 
one recovers the linear matrix equation for the scalar perturbation 



A 1 C , (scalar perturbation) 







d 2 W 

d(\nn) 2 



A 



dA- 



dW 



1 



din ft Shift H 2 fl 2 



D 2 W . (scalar perturbation) 



(4.17a) 



(4.17b) 



The case of a single scalar field has been much studied. Mukhanov et al |2(| derived the scalar perturbation 
equation by considering quadratic perturbations of the Einstein action. Hawking, Laflamme and Lyons p7{ utilized 
this equation in their discussion of the arrow of time. Using the Einstein field equations, Hwang |2£| gave an elegant 
interpretation in terms of uniform curvature gauge. Salopek and Stewart used the scalar equation to discuss 
microwave anisotropies and galaxy correlations arising from inflation. Dcruclle et al ]29| ] employed an equation which 
is basically equivalent to eq.( 4.17b ) in their discussion of extended inflation jl2| which incorporates two scalar fields. 
A path integral appro ach for multiple fields was considered in ref. |30|. However, the present form of the scalar 
perturbation equation ( 4.17b| ) possesses some distinct advantages over earlier formulations. 

Because there are no mass terms appearing on the right-hand-side, eq.(4.17b) is a particularly elegant form for the 
scalar perturbation equation. For example, in the large wavelength limit, one may neglect the last term appearing on 
the right-hand-side of eq.(4.17b). For an inflationary epoch, W then approaches a constant, 



^(ft,^^ 2 ) -> T/F°°(e,£> 2 ) 



(4.18) 



(i.e., it is independent of the integration parameter In ft). The numerical value of W°° and how it depends on D 2 are 
of high importance for observational cosmology. Approximate results will be given in the next section. 
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V. INITIAL CONDITIONS 



It is difficult to obtain analytic solutions of the scalar perturbation equation (4.17b) for two scalar fields. 



j perturbation equation ( jLiYop tor two scalar neids. Numerical 
computations have been considered in an earlier work |1C| . An approximate analytic method which is well known for 
the case of a single scalar field case will be extended to include multiple fields. Long wavelength and short w avelength 
perturbations will be matched at the Hubble radius. (Solutions of the tensor evolution equation (4.16b) and their 
applications to cosmology have already been examined, 010? an d they will not be discussed further in this section.) 

During a period of cosmological inflation, the initial conditions for the tensor perturbation y and the scalar per- 
turbation matrix W are set by quantum conditions: for wavelengths that are much shorter than the Hubble radius, 
1/H, one assumes that each mode began in the ground state. 

It is useful to rewrite the scalar evolution equation as a set of first order differential equations. One first introduces 
a momentum variable P, which is a matrix defined by 



P = A 



_! dW 



dlnQ 



The evolution equation for the pair (W, P) is then 

dW 



dP 



= AP, 
1 



D 2 A~ 1 W. 



(5.1) 

(5.2a) 
(5.2b) 



In general, (W, P) may be complex. 



A. Symmetry Conditions 



From the onset, it has been assumed that the scalar matrix S appearing in eq. (|4.15a| ) is symmetric. This requirement 
will impose several conditions on the fields appearing in the Riccati equation ( [4.17a]) . Since 



afb 2 = b 1 a2 , 

one may show that A^ 1 C is indeed symmetric. PW^ 1 is symmetric provided that 



W T P - P T W = . 



(5.3) 



(5.4) 



Using the first order evolution equation (5.2b), one ve rifie s that the derivative of the left-hand-side with respect to 
lnfi indeed vanishes. Hence, one typically imposes eq^ES^) as an initial condition. 



B. Conjugate Conditions 



In addition, it is possible to arrange that (W, P) are canonically conjugate. The conjugate conditions are 

WW^ -W*W T = 0, 

ppt _ p*pT = o ( 



WP f - W*P 7 



-il. 



(5.5a) 
(5.5b) 
(5.5c) 



where f represents the Hermitean conjugate, W' — (W*) T , and I is the 2-by-2 identity matrix, 7jj = Sij. If these 
conditions are met at one time, then they arc true for all times. For example, one may thus verify that the derivative 
of the left-hand-side with respect to In f2 vanishes (weakly) . 

Quantum considerati ons i mply the conjugate conditions, which legislate that (W, P) are intrinsically complex since 
i = \/^T appears in eq.( 5.5c ). The conjugate conditions will be used to prepare the cosmological system in the ground 
state at short wavelengths (se e Sec . C). 

The probability functional (2.3t) for scalar perturbations is determined to be 



V = exp - i J (fxf 1 ' 2 [D- 2 R, e] (WW f ) 1 [£>~ 2 i?, e\ J 



(5.6) 
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This functional is peaked about a flat and homogeneous field configuration: R = = D 2 e. For small deviations about 
this configuration, the 2-point correlation functions are found to be: 



< 



D- 2 R(x)] \D- 2 R(y)] >= /"^(x) [WW\x)] n 5 3 (x-y), 



< 



~D- 2 R(x)} e(y) >= r 1/2 (x) [WW\x)} u S 3 (x - y) , 
< e(x) e(y) >= f-^ 2 (x) [WW\x)] 22 S 3 (x - y) . 
(Eq.(pij|) follows from the Riccati equation (4.17a) and the conjugate conditions ( |5.5a| -c): 

1 



S-S* 



2D 4 



2D 4 



-A~ 



dW :W-l-^(W->) 



dlnQ 



pw^ (w^) 1 w- 1 - p*w t (w T ) 1 (w~ 1 ] 



etc.: 



in addition, A 1 C was assumed to be real. 



(5.7) 
(5.8) 
(5.9) 

(5.10) 
(5.11) 



C. Short- Wavelength Behavior 

One of the main assumptions of the inflationary model is that short wavelength fluctuations originated in the ground 
state. Analytic expressions are given for two scalar fields in such a state. 

The Laplacian operator D 2 for the conformal 3-metric fij appears throughout this work. Much of the subsequent 
operator analysis is simplified by introducing cigcnfunctions u(x,k) with eigenvalues — k 2 : 



D 2 u(x, k) = —k 2 u(x, k) . 

In the basis of the eigenvectors, the operator W describing the ground state |10| evolves according to 



W(n,e,k) = 



1 



Akr v — 1 



2kn 



X-\ k/(HCl)»l 



for wavelengths much shorter than the Hubble radius, k/(Hfl) >> 1. Conformal time r is given by 



dnn 2 



i 



H(fl,e) 



and the components of the transformation matrix X are 

Y 1 
A H = 7" 



4<91nfi 



Xl2 — 



de 



(5.12) 



(5.13) 



(5.14) 



(5.15a) 



X21 



1 



Xo 



4<9mft ! de 
X is actually the Jacobi matrix for the change of variables, (41nO,e) 



!>i, ^2 )■ The inverse X 1 



(x- 



= 4- 



(x- 



= 4- 



.d\iin 



(5.15b) 



(5.16a) 



{X ) 2 1 = -T— , (X ) 2 2 = -5— , 

0(pl 002 

is the Jacobi matrix for the inverse transformation, (0i, 02 ) — > (41nfi, e). Since 

1 



A 



1( v-l\T 



x- L (x 



(5.16b) 



(5.17) 



follows from eq.(4.15b), one may readily verify the symmetry condition (5.4) and the conjugate conditions ( |5.5a| -c). 
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D. Matching Long and Short Wavelength Fluctuations 



The short- wavelength approximation breaks down when the comoving wavelength equals the Hubble radius, 
k/{HVl) ~ 1, at which time W(VL, e, k) quickly approaches a constant value W°°(e,k) = W(oo,e,k) given approxi- 
mately by 



W°°(e,k) 



H 

TIP 



x- 1 



k=HQ, 



(5.18) 



W°°(e, k) contains virtually all of the observable properties of the cosmological model. 

For a single scalar field, <f>\ = <fi, one obtains the well-known expression for the variable zeta, Q: 



4 V2P d<t> 



k=iin 



3H 2 
V2P 



N 



k=HQ, 



The power spectrum for £ is 



k=HCt 



For two scalar fields, the power spectrum for Q is 



(5.19) 



(5.20) 



(5.21) 



Examples for two scalar fields will be given in the next section. 

For scalar perturbations, the resulting cosmic microwave background temperature anisotropy is related to zeta 
through 



AT(x) 
T 



-C(x)/15, |x| = 11, 700Mpc 



(5.22) 



where \x\ — 11, 700 Mpc is the comoving radius for the surface of last scattering, assuming that the present Hubble 
parameter is Hq = 50 km / s/ Mpc. 



VI. DOUBLE INFLATION POWER SPECTRUM 



As an example of the quadratic curvature approximation, one computes fluctuations arising from inflation that give 
rise to galaxy clustering as well as microwave background anisotropies. Consider the potential for two scalar fields, 



V(0 O ) = Wi) + V 2 (0a) 



(6.1) 



which is the sum of two separable potentials, Vi(cf>i) and V^ifo)- This potential has been used to describe the double 
inflation model, where <j>\ dominates the energy density of the first epoch of inflation, and (j>2 dominates the second 
epoch. 

Although exact solutions of the Hubble function H(cj) a ) exist |^], for illustrative purposes, I will be content to utilize 
the well-known slow-roll approximation, 



Vi(0i) + y 2 (02) 



1/2 



(6.2) 



where one neglects partial derivatives of the Hubble function in the SHJE ( ^.3b| ). This is one of those case s wer e 
finding a solution depending on two homogeneous parameters is very difficult, and I will instead integrate eq.(B.ll). 
One may verify the solution: 



In ft 



Jo 



y l dVi 

d<f,[ 



J 2 dV 2 ' 



dVi 



dV 2 



(6.3a) 
(6.3b) 
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The inverse transformation matrix is 



(x- 



: dVi 



dV 2 



(6.4a) 



(X x ) 2 i - g^" , (X : ) 

O0i 



1 

dV 2 



(6.5) 



In the slow roll approximation, k — Hfl may be approximated by k ex. Q. The power spectrum, V^(k), eq.(5.21), is 
then given by: 



4tt 2 



(Vi + v 2 ) 



Vx. 



tdVx 



dVi 



In [fc/tlO^Mpc" 1 )] = 60- 



01 ^Wi) 



dV 2 
(>d,' 



a</>' 



(6.6) 
(6.7) 
(6.8) 



It is assumed that the scale ko = 10 _4 Mpc experienced 60 e-foldings of inflation after it crossed the Hubble radius dur- 
ing inflation. For simplicity, e is assumed to be a homogeneous parameter. (Inhomogeneities in e lead to isocurvature 
perturbations which have been considered by Sasaki and Yokoyama [5jJ.) 

These analytic results verify the numerical simulations of ref . fl0|| . Heuristic arguments for these results have been 
given previously in refs. p2| and || . The double inflation model is strongly constrained by cosmological observations 



The case of quadratic scalar field potentials, 



Vv 







1 



2 ±2 

m\<p x , 



v 2 (<t> 2 ) 



- 2 j.2 



(6.9) 



is particularly easy to implement. In Fig.(l), power spectra for zeta are plotted for two models: inflation with a single 
scalar field, and inflation with two scalar fields. Both models have been normalized using COBE's 2-year data set [p4[ : 
CT s fcy(10 ) = [< (AT) 2 >] 1 / 2 = 30.5±2.7/iif. Using the methods of refs. 0], M, a small correction has been made for 
gravitational waves. For the single field model, the potential is V{4>) — m 2 (jr /2, with m = 5.45 x 10~ 6 (recall that 
mp = \Z8tt). For the double field model, the potential is given by eq.(|6j|), with mi = 1.435 x 10~ 5 and m 2 = mi/6. 

In Fig. (2), the observed galaxy data are displayed using the power spectrum for the density field, S = Sp/p. The 
data were compiled using several galaxy surveys p5| , and it has been assumed that infrared galaxies trace the mass 
distribution, i.e., the bias factor bjnAS = 1 is unity. For a single field, the theoretical curve is unsatisfactory because 
it substantially overestimates the power at short scales. For the double inflation model, the fit is satisfactory provided 
the parameters have been carefully chosen. Peter et al |33| have considered additional cosmological tests. These 
models are the subject of a continuing investigation. 



VII. CONCLUSIONS 



In elementary applications to mechanics, Hamilton- Jacobi theory is useful because it enables one to compute the 
constants of integration of a dynamical system. This useful feature becomes even more attractive when one considers 
Hamilton- Jacobi theory for general relativity. Einstein's field equations contain two constraints which legislate that 
the energy density and the momentum density vanish. These constraints are, in fact, constants of integration for the 
remaining evolution equations. Hence by solving the Hamilton- Jacobi equation for general relati vity, one solves for 
the constraints as well as the other constants of integration which include, for example, zeta, eq.( |5.21 ), for multiple 
scalar fields. 

Another attractive feature of HJ theory for general relativity is that it yields a covariant formulation of gravity. One 
may perform calculations without selecting an explicit choice for the temporal and spatial coordinates. In fact, the 
generating functional is manifestly invariant under transformations of the spatial coordinates since it is constructed 
by integrating the Ricci tensor and spatial derivatives of the various scalar fields over the 3-volume. The invariance 
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of the generating functional under a change of time hypersurface is more subtle. Using a spatial gradient expansion, 
one may express the generating functional as a line integral in superspace. Each contour of integration in the line 
integral corresponds to a specific choice of time hypersurface, but all choices yield the same answer. An integrability 
condition is met because the generating functional is invariant under spatial coordinate transformations. 

For two scalar fields, one employs the method of characteristics to transform the old variables (fa, fa, fij) to new 
ones (fi, e, fij). The integration parameter is arbitrary, but the choice of the conformal factor = fl(<p a ) simplifies the 
subsequent integration. The remaining new variables e and fij are constants in the long-wavelength limit. One may 
sum an infinite subset of terms in the spatial gradient expansion by making an ansatz for the generating functional 
which is quadratic in the fields D 2 e, R and Rij. 

HJ methods have proved particularly useful in applications to cosmology. As an illustration of the theory, the 
fluctuation spectra for double inflation was computed analytically confirming previous work. One may adjust the 
double inflation model to give an adequate explanation of galaxy clustering as well as large angle microwave background 
fluctuations. 



ACKNOWLEDGMENTS 



The author thanks J.M. Stewart and D. Page for useful discussions. This work was supported by the Natural 
Sciences and Engineering Research Council of Canada, and the Canadian Institute for Theoretical Astrophysics. 



[1] D.S. Salopek and J.M. Stewart, Phys. Rev. D 51, 517 (1995). 

[2] J. Parry, D.S. Salopek and J.M. Stewart, Phys. Rev. D 49, 2872 (1994). 

[3] D.S. Salopek and J.M. Stewart, Phys. Rev. D 47, 3235 (1993). 

[4] D.S. Salopek and J.M. Stewart, Class. Quantum Grav. 9, 1943 (1992). 

[5] D.S. Salopek, 43, 3214 (1991); ibid 45, 1139 (1992). 

[6] D.S. Salopek and J.R. Bond, Phys. Rev. D 42, 3936 (1990); ibid, 43, 1005 (1991). 

[7] K.V. Kuchaf, in Proc. of 4th Canadian Conference on GR and Relativistic Astrophys., May 16-18 1991, ed. G. Kunstatter 

et al, p. 211 (World Scientific, 1992). 
[8] G.F.R. Ellis, J.C. Hwang and M. Bruni, Phys. Rev. D 40, 1819 (1989). 

[9] D.S. Salopek, Phys. Rev. Lett, 69, 3602 (1992); Proc. of the International School of Astrophysics, "D. Chalonge", Third 
Course: Current Topics in Astrofundamental Physics, Erice, Italy, September 4-16, 1994, ed. N. Sanchez (World Scientific 
Publishers, 1995). 

[10] D.S. Salopek, J.R. Bond and J.M. Bardeen, Phys. Rev. D 40, 1753 (1989). 

[11] F. Lucchin and S. Matarrese, Phys. Rev. D 32, 1316 (1985); J. Yokoyama and K. Maeda, Phys. Lett. B207, 31 (1988). 
[12] D. La and P.J. Steinhardt, Phys. Rev. Lett. 62, 376 (1989); EW. Kolb, D.S. Salopek and M.S. Turner, Phys. Rev. D 42, 

3925 (1990). 
[13] Y.B. Zel'dovich, A & A, 5, 84 (1970). 

[14] K.M. Croudace, J. Parry, D.S. Salopek and J.M. Stewart, Ap. J. 423, 22 (1994). 

[15] D.S. Salopek, J.M. Stewart and K.M. Croudace, Mori. Not. Roy. Astr. Soc. 271, 1005 (1994). 

[16] S. Matarrese, O. Pantano and D. Saez, Phys. Rev. D 47, 1311 (1993). 

[17] E.J. Copeland, E.W. Kolb, A.R. Liddle and J.E. Lidsey, Phys. Rev. D 48, 2529 (1993). 

[18] L. Moscardini, S. Borgani, P. Coles, F. Lucchin, S. Matarrese, A. Messina, M. Plionis, Astrophys. J. 413, 55 (1993). 
[19] C. Misner, K.S. Thorne and J.A. Wheeler, Gravitation (New York: Freeman, 1973). 
[20] D.S. Salopek, Phys. Rev. D 46, 4373 (1992). 

[21] A.O. Barvinsky, Physics Reports 230, 5 & 6, 237 (1993); Phys. Lett. B241, 201 (1990); A.O. Barvinsky and A.Y. Kamen- 

shchik, Phys. Lett. B332, 270 (1994). 
[22] A. Peres, Nuovo Cim. 26, 53, (1962). 

[23] S. Tomonaga, in Quantum Electrodynamics, ed. J. Schwinger (Dover, New York, 1958); J. Schwinger, ibid.. 

[24] V. Moncrief and C. Teitelboim, Phys. Rev. D 6, 966 (1972). 

[25] D.S. Salopek, J.M. Stewart and J. Parry, Phys. Rev. D 48, 719 (1993). 

[26] V.F. Mukhanov, Zh. Eksp. Tor. Fiz. 94, 1 (1988) [Sov. Phys. JETP 67, 1297 (1988)]. 

[27] S. Hawking, R. Laflamme and G. Lyons, Phys. Rev. D 47, 5342 (1993). 

[28] J. Hwang, Phys. Rev. D 48, 3544 (1993). 

[29] N. Deruelle, C. Gundlach and D. Langlois, Phys. Rev. D 45, 3301 (1992). 



1G 



[30] S. Anderegg and V. Mukhanov, HEP-TH-9403091 preprint (1994). 
[31] M. Sasaki and J. Yokoyama, Phys. Rev. D 44, 970 (1991). 

[32] A. A. Starobinsky, Pis'ma Zh. Eksp. Tcor. Fiz. 42, 124 (1985) (JETP Lett. 42, 152 (1985)). 
[33] P.Peter, D. Polarski and A. A. Starobinsky, Phys. Rev. D 50, 4827 (1994). 
[34] C.L. Bennett et al, Astrophys. J., 436, 423 (1994). 

[35] J.A. Peacock and S.J. Dodds, Mon. Not. Roy. Astron. Soc. 267, 1020 (1994). 



17 



VIII. FIGURE CAPTIONS 



Fig.(l): Power spectra for zeta are shown for two models: (1) standard single field inflation and (2) double inflation. 
Both utilize quadratic scalar field potentials. The spectra are normalized using large angle microwave background 
measurements. A small correction has been made for gravitation radiation. Double inflation gives less power at short 
scales. 

Fig. (2): The power spectrum for the observed density field 5 = Sp/p is shown. Standard single field inflation gives 
a poor fit of the data at shorter distances, k < 10~ 1,4 Mpc -1 . By adjusting several free parameters, double inflation 
removes this problem. 
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POWER SPECTRA FOR ZETA 
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POWER SPECTRA FOR DENSITY PERTURBATION 
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